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Anisotropic flow of charged hadrons, pions and (anti-)protons measured 
at high transverse momentum in Pb-Pb collisions at v /s^ = 2.76 TeV 



The elliptic, V2, triangular, V3, and quadrangular, V4, azimuthal anisotropic flow coefficients are mea- 
sured for unidentified charged particles, pions and (anti-)protons in Pb-Pb collisions at <Js m = 



2.76 TeV with the ALICE detector at the Large Hadron Collider. Results obtained with the event 
plane and four-particle cumulant methods are reported for the pseudo-rapidity range |tj| < 0.8 at 
different collision centralities and as a function of transverse momentum, pj, out to pj = 20 GeV/c. 
The observed non-zero elliptic and triangular flow depends only weakly on transverse momentum for 
pj > 8 GeV/c. The small pj dependence of the difference between elliptic flow results obtained from 
the event plane and four-particle cumulant methods suggests a common origin of flow fluctuations 
up to pj — 8 GeV/c. The magnitude of the (anti-)proton elliptic and triangular flow is larger than that 
of pions out to at least pj = 8 GeV/c indicating that the particle type dependence persists out to high 
Pi- 
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The goal of ultra-relativistic nucleus-nucleus collisions is to study nuclear matter under extreme con- 
ditions. For non-central collisions, in the plane perpendicular to the beam direction, the geometrical 
overlap region, where the highly Lorentz contracted nuclei intersect and where the initial interactions 
occur, is azimuthally anisotropic. This initial spatial asymmetry is converted via interactions into an 
anisotropy in momentum space, a phenomenon referred to as transverse anisotropic flow (for a review 
see HI). Anisotropic flow has become a key observable for the characterization of the properties and the 
evolution of the system created in a nucleus-nucleus collision. 

Identified particle anisotropic flow provides valuable information on the particle production mechanism 
in different transverse momentum, pi, regions ID. For pj < 2 — 3 GeV/c, the flow pattern of different 
particle species is qualitatively described by hydrodynamic model calculations Q. At intermediate pj, 
3 < pj < 6 GeV/c, the observed flow of the baryons is larger than that of the mesons 13JIU. For pj > 8 
GeV/c, the fragmentation of high-energy partons, resulting from initial hard scatterings, is expected to 
play the dominant role. While traversing the hot and dense matter these partons experience collisional 
and radiative energy loss J5J|6], which are strongly dependent on the thickness of the created medium Q. 
In the azimuthally asymmetric system, the energy loss depends on the azimuthal emission angle of the 
partem, which leads to an azimuthal anisotropy in particle production at high px EEl. 

The magnitude of the anisotropic flow is characterized by the coefficients in the Fourier expansion of the 
azimuthal distribution of particles with respect to the collision symmetry plane ifTTjl fTTTl : 

v n (pT,n) = {cos[n(<l>-*i> n )}}, (1) 

where pj, t], and are the particle's transverse momentum, pseudo-rapidity, and the azimuthal angle, 
respectively, and *P n is the «-th harmonic symmetry plane angle. For a smooth matter distribution in the 
colliding nuclei, the symmetry planes of all harmonics coincide with the reaction plane defined by the 
beam direction and the impact parameter, the vector connecting the centers of the two colliding nuclei at 
closest approach. In this case, for particles produced at midrapidity, all odd Fourier coefficients are zero 
by symmetry. Due to event-by-event fluctuations of the positions of the participating nucleons inside 
the nuclei, the shape of the initial energy density of the heavy-ion collision in general is not symmetric 
with respect to the reaction plane, and the *P n may deviate from the reaction plane. This gives rise to 
non-zero odd harmonic coefficients lTT2l fT3l [141 [T31 [T6l [171 [T8ll . and contributes to the difference in flow 
coefficients calculated from two- or multi-particle azimuthal correlations, and also to the difference in v n 
measured with respect to different harmonic symmetry planes. 

Large elliptic flow, V2, and significant triangular flow, V3, were observed at the Relativistic Heavy Ion 
Collider (RHIC) Eg 123 ED and at the Large Hadron Collider (LHC) 1221125112311251125112711281. In this 
paper we present the measurement of unidentified charged particle anisotropic flow out to pj = 20 GeV/c, 
and for protons and charged pions Q out to pj = 16 GeV/c. We also present unidentified charged particle 
quadrangular flow, V4, measured with respect to the second OP2) and fourth (^4) harmonic symmetry 
planes. 

The data sample recorded by ALICE during the 2010 heavy-ion run at the LHC is used for this analysis. 
Detailed descriptions of the ALICE detector can be found in 112911301 |3~T1 . The Time Projection Chamber 
(TPC) was used to reconstruct charged particle tracks and measure their momenta with full azimuthal 
coverage in the pseudo-rapidity range |tj | < 0.8, and for particle identification via the specific ionization 
energy loss, dE/dx, in the transverse momentum region pj > 3 GeV/c. Two scintillator arrays (VZERO) 
which cover the pseudo-rapidity ranges —3.7 < r\ < — 1.7 and 2.8 < v\ < 5.1 were used for triggering, 
and the determination of centrality [32] and symmetry planes. The trigger conditions and the event 
selection criteria are identical to those described in Il22ll23l[32l . Approximately 10 7 minimum-bias Pb- 
Pb events with a reconstructed primary vertex within ±10 cm from the nominal interaction point in the 
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Fig. 1: (color online) V2, V3, and V4 measured for unidentified charged particles as a function of transverse mo- 
mentum for various centrality classes. The dashed line represents the WHDG model calculations for neutral pions 
V2 B31 extrapolated to the LHC collision energy. For clarity, the markers for V3 and V4/vp 2 results are slightly 
shifted along the horizontal axis. Error bars (shaded boxes) represent the statistical (systematic) uncertainties. 



beam direction are used for this analysis. Charged particles reconstructed in the TPC in \t]\ < 0.8 and 
0.2 < pt < 20 GeV/c were selected. The charged track quality cuts described in ll22l were applied to 
minimize contamination from secondary charged particles and fake tracks. The charged particle track 
reconstruction efficiency and contamination were estimated from HIJING Monte Carlo simulations ll33l 
combined with a GEANT3 [34] detector model, and found to be independent of the collision centrality. 
The reconstruction efficiency increases from 70% to 80% for particles with 0.2 < pj < 1 GeV/c and 
remains constant at 80 ± 5% for pj > 1 GeV/c. The estimated contamination by secondary charged 
particles from weak decays and photon conversions is less than 6% at pj = 0.2 GeV/c and falls below 
1% for p T > 1 GeV/c. 

The selection of pions and protons at pj > 3 GeV/c is based on the measurement of the dE/dx in the 
TPC, following the procedure described in Il3"5l . Enriched pi on (proton) samples are obtained by selecting 
tracks from the upper (lower) part of the expected pion (proton) dE/dx distribution. For example, protons 
were typically selected, depending on their momentum, in the range from to —3a or from —1.5a to 
—4.5a around their nominal value in dE/dx, where a is the energy loss resolution. Note that dE/dx 
of pions is larger than that of protons in the pj range used for this study. The track selection criteria 
have been adjusted to keep the contamination by other particle species below 1% for pions and below 
15% for protons. The pion and proton vi and V3 are not corrected for this contamination. The systematic 
uncertainties in V2 and V3 related to the purity of the pion and proton samples are 2% for pj < 8 GeV/c 
and 10% for p T >8 GeV/c. 

The flow coefficients v„ are measured using the event plane method (v„{EP} |Q]]) and the four-particle 
cumulant technique (v„{4} [363), which have different sensitivity to flow fluctuations and correlations 
unrelated to the azimuthal asymmetry in the initial geometry ("non-flow"). The non-flow contribution to 
v„{4} is estimated to be negligible from analytic calculations and Monte Carlo simulations ll3^ 13811391 . 
The contribution from flow fluctuations was shown to be negative for v„{4} and positive for v„{EP} (TJ. 

The orientation of the symmetry planes *P n is reconstructed from the azimuthal distribution of hits mea- 
sured by the VZERO scintillators. The large gap in pseudo-rapidity between the charged particles in 
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Fig. 2: (color online) Comparison of the ALICE results on v n (pj) obtained with the event plane method to the anal- 
ogous measurements from ATLAS 0261 and CMS [27 1 collaborations, as well as v? measurements by STAR 1441 . 
Only statistical errors are shown. 



the TPC and those in the VZERO detectors greatly suppresses non-flow contributions to the measured 
v„{EP}. Assuming that there is no anisotropic flow in pp collisions, the non-flow contributions can be 
estimated by comparing the azimuthal correlations measured in heavy-ion collisions to those in pp. It 
was observed that the two-particle azimuthal correlations in pp and the most peripheral Au-Au collisions 
at -y/^NN = 0.2 TeV are very similar BUI , which suggests that non-flow dominates correlations in the cen- 
trality range 80-90%. The systematic uncertainty from the remaining non-flow, 5,^ ent , in the measured 
v„{EP} coefficients was estimated based on the equation: 



scent _ 80-90% , / ™ n ^ 

°n ~ V n y MC ent ' ^ 

where v™- 90% and M 8 °- 90% are the magnitude of v„ and average multiplicity for the centrality range 80- 
90%, respectively, and M cent is the average multiplicity in a given centrality class. The non-flow increases 
with pj and from central to peripheral collisions. For example, the non-flow contributions to V2 in 5-10% 
(40-50%) most central collisions are about 1% (2%) at pj = 1 GeV/c and reach up to 10% (12%) for 
Pt > 10 GeV/c. Other sources of systematic uncertainties were evaluated from the variation of the results 
with different cuts on the reconstructed collision vertex and the centrality estimated from the charged 
particle multiplicity measured in the TPC and VZERO detectors. Changes due to variations of the track 
selection criteria and the difference of the results obtained using only positively or negatively charged 
particles were considered as a part of the systematic error. The difference in the extracted coefficients 
using one or the other of the two VZERO detectors was found to be below 1% for V2 and V3, and below 
5% for V4 over the measured region of transverse momentum. The combined results from correlations 
with both VZERO detectors are denoted as v„{EP, \Ar\ \ > 2.0} in the following. The contributions from 
all sources were added in quadrature as an estimate of the total systematic uncertainty. The resulting 
systematic uncertainties in V2 are 3% for 0.9 < pj < 1 GeV/c and (^2%) for 9 < pj < 10 GeV/c in 
the 5-10% (40-50%) centrality class. The resulting systematic uncertainties in V3 are 3% for 0.9 < pr < 1 
GeV/c and increase to 6% (10%) for 7 < p T < 9 GeV/c for centrality 5-10% (40-50%). We assign an 
8% (16%) systematic uncertainty to V4 for 0.9 < pj < 1 GeV/c in the 5-10% (40-50%) centrality class, 
while for pj > 6 GeV/c the systematic uncertainty is dominated by non-flow contributions. 
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Fig. 3: (color online) Relative event-by-event elliptic flow fluctuations for unidentified charged particles versus 
transverse momentum for different centrality classes. For clarity, the markers for centrality classes > 10% are 
slightly shifted along the horizontal axis. Error bars (shaded boxes) represent the statistical (systematic) uncertain- 
ties. 



Figure [T] shows unidentified charged particle vi, V3, and V4 as a function of transverse momentum for 
different centrality classes. The difference between V2{EP} and v%{\} for pj < 7 GeV/c is predominantly 
due to flow fluctuations. The measured V2 at pj > 8 GeV/c is non-zero, positive and approximately 
constant, while its value increases from central to mid-peripheral collisions. The observed V2{EP} at 
pj > 10 GeV/c is fairly well described by extrapolation to the LHC energy iHTI of the WHDG model 
calculations Il42l for V2 of neutral pions including collisional and radiative energy loss of partons in 
a Bjorken-expanding medium P31 . The coefficient V3 exhibits a weak centrality dependence with a 
magnitude significantly smaller than that of v%, except for the most central collisions. Unlike V3, which 
originates entirely from fluctuations of the initial geometry of the system, V4 has two contributions, which 
are probed by correlations with the ^2 and ^4 symmetry planes. The measured V4/>j< 4 {EP} does not 
depend strongly on the collision centrality which points to a strong contribution from flow fluctuations. 
In contrast, v 4 /^> 2 {EP} shows a strong centrality dependence which is typical for correlations with respect 
to the true reaction plane. The difference between the two, indicative of flow fluctuations, persists at least 
up to pi = 8 GeV/c. 

Figure [2] compares our results obtained with the event plane method for 30-40% centrality to the anal- 
ogous measurements by ATLAS G61 and CMS ETl collaborations, and results obtained at RHIC by 
the STAR 043 collaboration. An excellent agreement is observed between results from all three LHC 
experiments. V2(pt) at top RHIC energy has a peak value about 10% lower than at LHC although is very 
similar in shape. 

To investigate further the role of flow fluctuations at different transverse momenta we study the relative 
difference between v 2 {EP} and v 2 {4}, [(v 2 {EP} 2 - v 2 {4} 2 )/(v 2 {EP} 2 + V2I4} 2 )] 1 / 2 , which for small 
non-flow is proportional to the relative flow fluctuations o V2 /{v2) CD- Figure [3] presents this quantity 
as a function of transverse momentum for various centrality classes. The relative flow fluctuations are 
minimal for mid-central collisions and become larger for peripheral and central collisions, similar to 
those observed at RHIC energies HI. It is remarkable that in the 5-30% centrality range, relative flow 
fluctuations are within errors independent of momentum up to pj ~ 8 GeV/c, far beyond the region where 
the flow magnitude is well described by hydrodynamic models (pj < 2 — 3 GeV/c). This indicates a 
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Fig. 4: (color online) Unidentified charged particle V2, V3, and V4 integrated over the transverse momentum range 
10 < pj < 20 GeV/c as a function of collision centrality, with the more central (peripheral) collisions shown 
on the left-(right-)hand side, respectively. The dashed line represents the WHDG model calculations for neutral 
pions [43 j extrapolated to the LHC collision energy. Error bars (shaded boxes) represent the statistical (systematic) 
uncertainties. 



common origin for flow fluctuations, which are usually associated with fluctuations of the initial collision 
geometry, at least up to the regime where hard scattering and jet energy loss are expected to dominate. 
The ratio develops a momentum dependence, starting to increase at pj ~ 1.5 GeV/c, for more peripheral 
collisions (30-50%), and in most central collisions (0-5%), where it is most pronounced. In both cases, 
the relative contribution of non-flow effects is expected to be the largest. 

Figure |3J shows unidentified charged particle vi, V3, and V4 averaged over 10 < pj < 20 GeV/c as a func- 
tion of centrality. V2 increases from central to peripheral collisions. No significant difference between 
V2{EP} and V2{4} results is observed, which might indicate that the fluctuations of the initial collision 
geometry become unimportant for pj > 10 GeV/c. The centrality dependence of V3 differs significantly 
from that of vi. V4 measured with respect to the second and fourth harmonic symmetry planes is con- 
sistent with zero within relatively large uncertainties. All these observations indicate that for pj > 10 
GeV/c the effect of fluctuations of the initial collision geometry might be very different compared to that 
at low and intermediate pj. 

Figure |5]pi'esents charged pion and proton v% and V3 as a function of pi in the 10-50% centrality range 
from the event plane method. The proton V2 and V3 are higher than that of pions out to pi = 8 GeV/c 
where the uncertainties become large. This behavior is qualitatively consistent with a picture where 
particle production in this intermediate pj region includes interaction of jet fragments with bulk matter, 
e.g. as in model BBl . The magnitude of the measured charged pion elliptic flow for pj > 8 GeV/c 
is compatible with that for unidentified charged particles, and 71° measured by PHENIX H6l in Au- 
Au collisions at y/s^N = 0.2 TeV, and reproduced by the WHDG model calculations for V2 of neutral 
pions 1351 . 

In summary, we have presented elliptic, triangular, and quadrangular flow coefficients measured by the 
ALICE collaboration in Pb-Pb collisions at y/s^ = 2.76 TeV over a broad range of transverse momen- 
tum. For pi > 8 GeV/c, we find that the unidentified charged particle vi and V3 are finite, positive and 
only weakly dependent on transverse momentum, while V4 is consistent with zero within rather large sta- 
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Fig. 5: (color online) vi (top) and V3 (bottom) of charged pion and proton as a function of transverse momentum 
for 10-50% centrality class compared to unidentified charged particles results from the event plane method. For 
clarity, the markers for v 2 and V3 at pj > 8 GeV/c are slightly shifted along the horizontal axis. PHENIX 7T° v 2 



measurements 1461 are also shown. The dashed line represents the WHDG model calculations for neutral pions | 
extrapolated to the LHC collision energy for the 20-50% centrality range. Error bars (shaded boxes) represent the 
statistical (systematic) uncertainties. 



tistical and systematic uncertainties. The observed difference in the centrality dependence of v^/\^ A and 
V4/»p 2 , and the results on V2 obtained with the event plane and four-particle cumulant methods indicate 
that the effect of flow fluctuations extends at least up to pi = 8 GeV/c and does not change significantly in 
magnitude. It shows that the effect of fluctuations of the initial collision geometry on particle production 
is similar at low and intermediate pj regions, which are considered to be dominated by hydrodynamical 
flow and quark coalescence, respectively. For pj > 10 GeV/c, where particle production is dominated 
by fragmentation of hard partons, the response to fluctuations of the initial collision geometry might be 
different, but more data is needed to study this regime in more detail. The pion V2 at LHC energies is very 
close to that measured at RHIC out to pi = 16 GeV/c and is reproduced by WHDG model calculations 
for pj > 8 GeV/c. The proton V2 and V3 are finite, positive, and have a larger magnitude than that of the 
pion for pj < 8 GeV/c, indicating that the particle type dependence, which is typical at low pj, persists 
out to high transverse momenta. 
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